LAMMPS (18 Feb 2013)
# Pour granular particles into chute container, then induce flow

atom_style	sphere
boundary	p p fm
newton		off
communicate	single vel yes

region		reg block -10 10 -10 10 -0.5 16 units box
create_box	1 reg
Created orthogonal box = (-10 -10 -0.5) to (10 10 16)
  1 by 1 by 1 MPI processor grid

neighbor	0.2 bin
neigh_modify	delay 0

# IMPORTANT NOTE: these values of stiffness (2000) and timestep (0.001)
# are used in this example file to produce a quick simulation and movie.
# More appropriate values for realistic simulations would be
# k = 2.0e5 and dt = 0.0001, as in bench/in.chute.

pair_style      gran/hooke/history 2000.0 NULL 50.0 NULL 0.5 0
pair_coeff	* *

timestep	0.001

fix		1 all nve/sphere
fix		2 all gravity 1.0 spherical 0.0 -180.0
fix		zlower all wall/gran 2000.0 NULL 50.0 NULL 0.5 0 		zplane 0.0 2000.0

region		slab block -9.0 9.0 -9.0 9.0 10.0 15.0 units box
fix		ins all pour 3000 1 300719 vol 0.13 50 region slab
Particle insertion: 402 every 3162 steps, 3000 by step 22135

compute		1 all erotate/sphere
thermo_style	custom step atoms ke c_1 vol
thermo		1000
thermo_modify	lost ignore norm no
compute_modify	thermo_temp dynamic yes

#dump		id all atom 1000 dump.pour

#dump		1 all image 1000 image.*.jpg type type #		axes yes 0.8 0.02 view 80 -30
#dump_modify	1 pad 5

run		25000
Memory usage per processor = 9.46169 Mbytes
Step Atoms KinEng 1 Volume 
       0        0           -0            0         6600 
    1000      402    779.98578            0         6600 
    2000      402    1424.4252            0         6600 
    3000      402    1425.0863      12.9927         6600 
    4000      804    1724.8532    47.325647         6600 
    5000      804    1598.6681    64.779741         6600 
    6000      804    1364.8231     61.75911         6600 
    7000     1206    1519.7287    69.028648         6600 
    8000     1206    1409.6162    62.604936         6600 
    9000     1206    1243.5299    53.012812         6600 
   10000     1608    1400.2632    57.300035         6600 
   11000     1608    1226.4327    53.185401         6600 
   12000     1608     1190.754    41.207269         6600 
   13000     2010    1309.6301    51.508727         6600 
   14000     2010      1182.74    46.857423         6600 
   15000     2010    985.17177    43.141534         6600 
   16000     2412    1041.1301    47.314349         6600 
   17000     2412    1038.3663    39.508439         6600 
   18000     2412    916.06571    39.312699         6600 
   19000     2814    1021.9901    42.073133         6600 
   20000     2814      959.466    41.811414         6600 
   21000     2814     761.2927    41.364383         6600 
   22000     2814    521.07124    40.192029         6600 
   23000     3000    423.49186     37.34074         6600 
   24000     3000    362.28297    22.586628         6600 
   25000     3000    223.45694    22.976343         6600 
Loop time of 9.28813 on 1 procs for 25000 steps with 3000 atoms

Pair  time (%) = 5.36181 (57.7276)
Neigh time (%) = 1.04765 (11.2794)
Comm  time (%) = 0.25077 (2.6999)
Outpt time (%) = 0.000460148 (0.00495415)
Other time (%) = 2.62743 (28.2881)

Nlocal:    3000 ave 3000 max 3000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    669 ave 669 max 669 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    13650 ave 13650 max 13650 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 13650
Ave neighs/atom = 4.55
Neighbor list builds = 1153
Dangerous builds = 0

unfix		ins
fix		2 all gravity 1.0 chute 26.0
run		25000
Memory usage per processor = 10.2619 Mbytes
Step Atoms KinEng 1 Volume 
   25000     3000    223.45694    22.976343         6600 
   26000     3000      113.473    17.666374         6600 
   27000     3000    102.76936    12.540567         6600 
   28000     3000    179.07068    11.501122         6600 
   29000     3000    289.11251    14.028808         6600 
   30000     3000    461.92265    17.681677         6600 
   31000     3000    710.67032    23.209353         6600 
   32000     3000    991.22981    30.179598         6600 
   33000     3000    1366.0205    36.321028         6600 
   34000     3000    1817.0366    45.963988         6600 
   35000     3000    2329.5052    51.957397         6600 
   36000     3000    2938.7862    59.654677         6600 
   37000     3000     3662.828    70.399596         6600 
   38000     3000    4375.2295    78.906063         6600 
   39000     3000    5250.4031    92.018452         6600 
   40000     3000    6032.6969    100.23325         6600 
   41000     3000    7053.8661    115.73204         6600 
   42000     3000    8063.1636    122.38643         6600 
   43000     3000    9181.8348    143.39138         6600 
   44000     3000     10270.31    149.34715         6600 
   45000     3000     11434.53    161.47217         6600 
   46000     3000    12731.277    171.19811         6600 
   47000     3000    13938.907    170.91903         6600 
   48000     3000    15244.449     192.6018         6600 
   49000     3000     16627.76     195.3851         6600 
   50000     3000    17863.962    213.64088         6600 
Loop time of 18.2656 on 1 procs for 25000 steps with 3000 atoms

Pair  time (%) = 11.864 (64.9527)
Neigh time (%) = 1.57208 (8.60681)
Comm  time (%) = 0.539768 (2.95511)
Outpt time (%) = 0.000703335 (0.0038506)
Other time (%) = 4.28904 (23.4816)

Nlocal:    3000 ave 3000 max 3000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    755 ave 755 max 755 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    13875 ave 13875 max 13875 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 13875
Ave neighs/atom = 4.625
Neighbor list builds = 803
Dangerous builds = 0
